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ABSTRACT 



We present the first analysis of the X-ray warm absorber and nuclear obscuration in the Seyfert 1.8 galaxy ESO 113-GOlO. We 
used archival data from a 100 ks XMM-Newton observation made in 2005. From high resolution spectroscopy analysis of the RGS 
data, we detect absorption lines originating from a warm absorber consisting of two distinct phases of ionisation, with logf « 3.2 
and 2.3 respectively. The higher-ionised component has a larger column density and outflow velocity {N^ x L6 x 10^- cm"-, v » 



-1 100 km s ') than the lower-ionised component (A'h ~ 0.5 x 10^^ cm 



-700 km s ). The shape of the optical-UV continuum 



and the large Balmer decrement (Ha/Hfj ~ 8) indicate significant amount of reddening is taking place in our line of sight in the host 
galaxy of the AGN; however, the X-ray spectrum is not absorbed by cold neutral gas intrinsic to the source. We discuss different 
explanations for this discrepancy between the reddening and the X-ray absorption, and suggest that the most likely solution is a dusty 
warm absorber. We show that dust can exist in the lower-ionised phase of the warm absorber, which causes the observed reddening 
of the optical-UV emission, whereas the X-rays remain unabsorbed due to lack of cold neutral gas in the ionised warm absorber. 
Furthermore, we have investigated the uncertainties in the construction of the Spectral Energy Distribution (SED) of this object due 
to obscuration of the nuclear source and the effects this has on the photoionisation modelling of the wam absorber. We show how 
the assumed SEDs influence the thermal stability of each phase and whether or not the two absorber phases in ESQ 113-GOlO can 
co-exist in pressure equilibrium. 

Key words. X-rays: galaxies - galaxies: active - galaxies: Seyfert - galaxies: individual: ESO 1 13-GOlO - techniques: spectroscopic 
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1 ESO 113-GOlO was first identified and catalogued as a galaxy 

■ (SBa) in the ESO/Uppsala Survey of the ESO(B) Atlas (.Lauberts, 
I Il982h . based on observations made with the ESO 1 m Schmidt 

telescope at La Silla, Chile. From its first-ever X-ray observa- 
tion with ROSAT in 1995 and follow-up optical spectroscopy in 
\ 1996 with the 2.2m ESO/MPG telescope at La Silla observatory, 

■ iPietsch et af] ([1998 ) have classified it as a Seyfert 1.8 galaxy at 
; a redshift of 0.025701. 

The next time ESO 113-GOlO was observed in the X-rays 
was with XMM-Newton in May 2001. From spectral analysis of 
! the 4 ks EPIC-pn data, iPorquet et aP (|2004 reported the pres- 
' ence of a soft X-ray excess and a highly redshifted Fe Ka line 
at 5.4 keV. The only other X-ray observation of ESO 1 13-GOlO 
to this date was 100 ks long and perforrn ed with XMM-Newton 
in November 2005. IPorquet et all (l2007l) reported strong rapid 
variably from their power spectral density timing analysis of the 
EPIC-pn data. Furthermore, from spectral analysis of the iron 
line band, no redshifted Fe Ka line was detected contrary to the 
2001 findings, while the presence of two narrow emission lines 
at about 6.5 keV and 7 keV was reported. 

In this work we present the first analysis of the RGS spec- 
tra from the 100 ks XMM-Newton observation of 2005, which 
show clear signs of warm absorber outflows. We also used the 
simultaneous EPIC-pn and OM data to aid us in the photoioni- 
sation modelling of the outflows. The structure of this paper is 
as follows. Section|2]describes the observations and data reduc- 
tion. Section |3] focuses on the nuclear obscuration of the AGN. 
Spectral analysis and photoionisation modelling of the warm ab- 



sorber are described in Sect. Construction of the SEDs and 
modelling with different SEDs are described in Sect. |5] We dis- 
cuss our findings in Sect.|6]and give concluding remarks in Sect. 

m 

Note that all the spectra and SEDs displayed in this paper are 
in the observed frame. The RGS spectra shown are background- 
subtracted. 



2. Observations and data reduction 

The ID of the analysed XMM-Newton observation is 
0301890101, stalling at 22:17 UTC on 10 November 2005 
with an EPIC-pn duration of about 102 ks. We extracted the 
RGS, EPIC-pn and OM data from the XMM-Newton Science 
Archive. All the observation data files were processed using 
SAS vlO.O. 



2.1. RGS data 

The RGS instruments were operated in the standard Spectro +Q 
mode. We filtered out time intervals with background count rates 
> 0.1 ct s ' in CCD number 9, which is the most affected by 
background flares. We removed about 8 ks, leaving a net ex- 
posure time of about 95 ks for each RGS instrument. The data 
were processed through the rgsproc pipeline task; the source 
and background spectra were extracted and the response matri- 
ces were generated. Using the rgscombine task, we combined 
the RGS 1 and RGS2 first-order spectra and response matrices in 
order to achieve better signal-to-noise for the purpose of spectral 
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fitting. The channels in the RGS spectrum were then re-binned 
by a factor of 3 to improve statistics. Finally, the modelling of 
the 6-38 A part of the RGS spectrum was performed using the 
SPEX v2.02.043 fitting package. 

2.2. EPIC-pn data 

During the XMM-Newton observation, the EPIC-pn camera was 
operated in the Full-Frame mode with the medium-filter applied. 
We extracted a single event, high energy light curve, with the 
#XMMEA_EP selection attribute, in order to create a set of Good 
Time Intervals (GTI) excluding intervals of flaring particle back- 
ground with count rates exceeding 0.4 ct s"'. We then applied 
the GTI to the production of the spectra, for which X-ray events 
with the FLAG -- attribute were selected. The spectra were 
extracted from a circular region of 35" radius centred on the 
source. The background was extracted from a nearby source- 
free region of the same size on the same CCD chip. The EPIC- 
pn data showed no evidence of pile-up, thus single and double 
events were selected. Response matrices were generated using 
the SAS tasks arfgen and rmfgen. The net exposure time of 
the EPIC-pn data used for spectral analysis was 77 ks. 

2.3. OM data 

The OM data from Image-mode operations were taken with U, 
UVWl, UVM2 and UVW2 filters. The OM data were processed 
with the omi chain pipeline. We performed aperture photometry 
on each image in a fully interactive way using the omsource 
program. The OM count rates were extracted from a circle of 12 
pixels radius (5.8") centred on the source nucleus. The source is 
point-like within the extraction region in the UV OM data. The 
background was extracted from a source-free region of the same 
radius. We then applied all the necessary corrections, i.e. for the 
point spread function (PSF) and coincidence losses, including 
time-dependent sensitivity (TDS) corrections. 



3. Intrinsic obscuration of tlie nuclear source 

3.1. Nuclear obscuration from the Balmer decrement 

ESO 113-GOlO is classified as a Seyfert 1.8 galaxy (with 
strong narrow components and very weak but still visible broad 
components of Hg /16563 and li/ 3 A4S61 emission lines) by 
iPietsch et al] (Il998h . following the IOsterbroc3 (11989 ') observa- 
tional classification of AGN. The Hg/H8 lin e ratio is not one of 
the properties given in iPietsch et al. 1(^1998), so we re-analysed 
their optical spectrum in order to obtain the Balmer decrement. 
We fitted both the Ha and H/3 lines with two narrow and broad 
Gaussian-profile components. The two forbidden lines of [N ii] 
/1/16548, 6583, which are blended with Ha, were also mod- 
elled with narrow Gaussian-profile components. To calculate the 
Ha/Hfi ratio we used the total flux of each line (including both 
narrow and broad components) since the broad component of the 
Hy8 is very weak and de-blending of the narrow and broad com- 
ponent of the lines increases the uncertainty in the flux measure- 
ments. So the Balmer decrement and consequently the reddening 
calculated here is an average over the NLR and BLR. We find 
the observed Balmer decrement Ha/Hj3 - 7.96. The FWHM 
is about 2000 km s"' for the broad components and about 500 
km s ' for the narrow components of the above lines. 
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To convert the Balmer decrement into reddening E(B - V) , 
we used the mean extinction curve of iGaskell & Benken (l2007h 
obtained for a sample of AGN using HST. Their best deter- 
mined extinction curves are flatter in the far-UV and are miss- 
ing t he /t2175 bump c o mpared to the stand ard Galactic curve 
(e.g. lOsterbrockl Il989l ICardelh et al.l Il989l) . The relationship 
/.i.obs = /i,in, 10-°-4^' gives 

log j = -0.4 (A,,H„ - A,,n8) (1) 

where Robs and Ri„t are the observed and intrinsic Balmer decre- 
ments Ho'/HyS; A^ Ha and A^i h/j represent the extinction at the 
wavelengths of the Balmer li nes. Using the parameteri sation of 
the mean extinction curve of iGaskell & Benkerl (l2007l) . we ob- 
tain the colour excess 

£(B-y) = 1.790 log (2) 

The scalar specifying the ratio of total to selective extinction 
Rv = Av/E(B - V) was fixed at 3.1. 

From theoretical recombination studies, the Balmer decre- 
ment Ha/HjS - 2.85 for Case B recombination at a tem- 

S ierature T = 10'^ K and electron density - 10'* cm"^ 
Osterbrock & Ferlandll2006 l). Case B refers to the recombina- 
tion which takes place in a typical nebula with a large enough 
optical depth that Lyman-line photons are scattered n times and 
are converted (if n > 3) into lower-series photons plus either Ly- 
a or two-c ontinuum photons, so that th ey cannot escape from 
the nebula (O sterbrock & Ferlandll2006 '). For the NLR and BLR 
of AGN, which have higher electron density than a nebula, this 
value requires modification due to contribution to the Ha line 
from collisional excitation. For the NLR of AGN the intrinsic 
Balmer decrement Ha/H/3 is about 3.1 (lOsterbrock & Ferlandl 
I2006h . This value however may be difl'erent for the BLR as 
clouds have higher density in this region. Ferguson & Ferlan'3 
(Il997h have computed the ratio of Balmer lines relative to H/3 for 
realistic conditions in BLR clouds. For a BLR with T - 10"* K 
and He = 10'^ cm ~^, they obtain an intrinsic Ha/Hj3 - 3.28. 
From observations, IWard et aP (1 19881) have measured the aver- 
age intrinsic Ha/H/3 in BLR to be 3.5, independent of any atomic 
physic s assumptions. Using a similar approach, ICarrera et al.l 
(12004 *) have also measured the intrinsic Ha/H/3 to be 3.43. 

So in our calculations we take Ri^t =3.1-3.5 and /?obs to 
be 7.96, which yields intrins ic E(B - V) = 0.64-0.73. We then 
use the extinction curve of iGaskell & Benker (2007) to correct 
all the optical-UV data for reddening in the host galaxy of the 
AGN, using E(B - V) = 0.64 and E(B - V) = 0.73 for two 
possible scenarios. 

3.2. Nuclear obscuration from the optlcal-UV continuum 

Intrinsic reddening will also affect the optical/UV continuum 
and change the shape of the so-called 'big-blue-bump', which 
was firs t attribu t ed to thermal emission from the accretion 
disc by I Shieldsl d 19781) . The optical/UV continuum is AGN- 
dominated in ESO 1 13-GOlO; for the purpose of our study con- 
taminating emission by the stars in the host galaxy would have a 
negligible effect on the SED and photoionisation modelling. We 
corrected the continuum in such a way that (1) the optical slope 
(ffopt) and (2) the relation between optical-to-X-ray spectral in- 
dex (ffox) and the logjo of the monochromatic optical luminosity 
at 2500 A (/opt) are consistent with AGN samples found in sur- 
veys. 
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lYoung et alJ (l2010l) have cross-correlated the Sloan Digital 
Sky Survey (SDSS) DR5 quasar catalog with the XMM-Newton 
archive, and have obtained aopt and aox for a sample of 327 
quasars with high X-ray signal-t o-noise ratio, where both optical 
and X-ray spectra are available. lYoung et all (1201 Ol) find ffopt - 
- 0.40, which is also simila r to the slope aopt = -0.44 found 
bv lVanden Berk et"al JipooT) for mean composite quasar spectra 
using a dataset of over 2200 spectra from the SDSS. For the un- 
corrected continuum, we find ffopt = -1.74, w hich indicates red- 
dening . So, we used the extinction curve of iGaskell & BenkeJ 
(l2007h to correct the optical/UV data for reddening to match the 
^opt - -0.40. We find that the amount of reddening required for 
this correction corresponds to E{B - V) - 0.39, which is less 
than E{B - V) range of 0.64 to 0.73, inferred from the Balmer 
decrement. 

Furthermore, we checked our values of anx an d Zopt against 
the aox-/opt relation found by lYoung etakl (l2010h . For the un- 
corrected case we find a ox - -1-2 and log /opt = 28.0 
(erg s"' Hz"'), and for the continuum corrected case with aopt - 
-0.40, we find ffox = -1.6 and log /opt = 29.1 (ergs"' Hz"'). In 
both cases, the values we find are consistent within the disper- 
sion in the Cox-^opt relation: for log /opt - 28.0, ffox is between 
-1.4 and -1.1; and for log /opt = 29.1, aox is between -1.6 and 
-1.3. On the other hand, for the Balmer-decrement corrected 
case, aox (-2.0 to -1.9) is not consistent with the aox of -1.7 to 
-1.3 from the aox-/opt relation for log /opt between 29.7 and 29.9. 
This may indicate that the Balmer decrement over-estimates the 
amount of reddening; we return to this issue in Sect. 16.11 



3.3. Intrinsic X-ray obscuration of the nuclear source by cold 
neutral gas? 

From the intrinsic reddening E{B - V) values calculated in Sect. 
13.11 from the Balmer decrement and in Sect. 13.21 based on the 
continuum reddening, the associated column density of cold 
gas in ESO 113-GOlO can be estimated. The relationship be- 
tween hydrogen column density A^h and optical extinction Ay 
can be written as A^h (cm"^) = a x 10^' Ay (mag), where the 
factor a is reported to be 2.22 i n Xj orenstein (1975), 1.79 in 
|Predehl&Sch mitt(1995), 1.89 in loTterbrock & Ferlandl(i2006b 
and 2.21 in lGiiver & Ozel(.2009;) . Therefore, one can calculate a 
range of values for A^h using the different values of E{B - V) and 



a given above. This gives A^h =(3.5-5.0) xlO cm"- if using 
E(B - V) calculated from the Balmer decrement and A^h = (2.2- 
2.7) xlO^' cm"^ if using E{B-V) calculated from the continuum 
reddening. 

However, analyses of the X-ray spectra show that ESO 113- 
GOI O is not intrinsically absorbed by a large column of neutral 
gas. lPorquet et al.l ( 12004 ^) obtain an upper limit of 8 x lO'^ cm"^ 
to intrinsic absorption from analysis of the 2001 EPIC-pn spec- 
trum. This is indeed much lower than the predicted values of 
(3.5-5.0) xlO^i cm"2 or (2.2 -2.7) xlO^' cm"^ cal culated above 
from the intrinsic reddening. iPorquet et al.l (l2007l) also mention 
that during the 2005 XMM-Newton observation, instead of in- 
trinsic neutral absorption expected below about 1 keV, a soft 
X-ray excess is observed. From analysis of the 2005 EPIC-pn 
spectrum, we find an upper limit of 8.9 x lO'** cm"^ (at 99% 
confidence level) to the A^h of intrinsi c neutral gas. This up per 
limit is compatible with that found by 'Porauet et alj (l2004l) for 
the 2001 data. We discuss this apparent discrepancy between 
the optical-UV reddening and X-ray absorption, and its impli- 
cations, in Sect. 16.1! 



4. Photoionisation modelling and spectral analysis 
of the warm absorber 

We began our modelling of the soft X-ray spectrum of ESO 113- 
GOIO by fitting the RGS continuum with a cosmologically red- 
shifted power-law. In all our fits the effects of the Galactic 
neutral absorption were included by applying the hot model 
(collisional ionis ation equihbrium ) in SPEX. This absorption 
model uses the IVerner et al.l (1199 6) cross-sections; for informa- 
tion about the atomic dat abase used in SPEX s ee the SPEX ref- 
erence manual. Assuming lLodders et al.l (l2009l) abundances, the 
Galactic Hi column density in our line of sight was fixed at 
A^H = 2.78 X 10-" cm"2 (Dickey & Lockman 1990) and the gas 
temperature at 0.5 eV to mimic a neutral gas. We note here that 
as requested by the referee of this paper, we also tested the tb- 
new model (Wilms et al. 201 1) for X-ray absorption in the ISM 
using the Wilms et al. (2000) abundances. We find the features 
in the vicinity of the K-edges of Ne and O and L-edge of Fe to be 
very similar in both the tbnew and hot models as shown in Fig.[T] 
Therefore, using the tbnew model would not change the results 
of our warm absorber study. Furthermore, we checked the im- 
pact of using the slightly smaller A^h value of 2.08 x 10"" cm"^ 
from the Leiden/ Arg entine/Bonn (LAB) Survey of Galactic H i 
(iKalberla et al.l2005l) on our analysis; we find that all parameters 
derived in this work remain unchanged within the given errors. 

The power-law fit, letting photon index and normalisation 
free, led to a reduced Chi-squared ixl) value of 1.6 for 1054 de- 
grees of freedom (d.o.f.). The fit was not satisfactory as there 
were absorption features, which were not fitted by the Galactic 
absorption model. At this stage we included one xabs compo- 
nent in the model. For an assumed SED of the source, xabs 
calculates the transmission through a slab of material where 
all ionic column densities are linked in a physically consistent 
fashion through a photoionisation model. The ionisation balance 
calculations were performed using version C08.00 of Cloud>0 
(Ferland et al. 1998) with Lodders et al. (2009) abundances. For 
an assumed SED, we made runs with Cloudy for a grid of ion- 
isation parameter (^) values (log^ between -8.5 and -1-6.5 with 
steps of 0.1) in order to calculate the equilibrium i on concen- 
trations. The ionisation parameter ^ (introduced by Tart er et aH 
1969), which is used to describe photoionisation equilibrium, is 
defined as 



L 



(3) 



where L is the luminosity of the ionising source over the 1- 
1000 Ryd band (in erg s"'), n the hydrogen number density (in 
cm"^) and r the distance between the ionised gas and the ionising 
source (in cm). Therefore, ^ is in units of erg cm s"' . 

Using the auxiliary program xabsinput in the SPEX pack- 
age, the Cloudy output is converted into an input file, containing 
the temperature and ionic column densities as a function of ^, for 
the xabs photoionised absorption model of SPEX. In the mod- 
elling we fitted the ionisation parameter (^), the equivalent hy- 
drogen column density (A^h) of the warm absorber, its flow and 
RMS velocities. As a start, we adopted the SED of the Seyfert-1 
galaxy Mrk 509 (but see Sect.|5]for detailed modelling with dif- 
ferent SEDs for ESO 1 13-GOlO). This SED was established with 
the help of simultaneous data from a large multi-wavelength 
campaign (Kaastra et al. 201 lb, and broad-b and modelling of the 
underling continuum (Mehdipou r et al.ir201 1), and is likely to be 
representative of the real SED shape for a typical Seyfert galaxy 
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Fig. 1. The tbnew and hot ISM X-ray absorption models as described in Sect. |4]applied to the power-law continuum of ESO 113- 
GOIO, shown in the vicinity of the K-edges of Ne (left panels) and O (right panels) and L-edge of Fe (middle panels). 



(SED G in Fig. |4]l. The inclusion of the xabs component \m- 
proved;^^^ to 1.22 (1050 d.o.f.). There were, however, absorption 
lines in the lower energy part of the spectrum which could not be 
fitted with a single xabs phase, since they are less highly ionised 
and have a different velocity. So we introduced a second xabs 
phase in the model to obtain a better fit to the data. The addition 
of the second phase further improved to 1.12 (1046 d.o.f.) and 
all the significantly detected absorption features were fitted. We 
tested adding another xa^Ji component, but this did not make the 
fit any better, so we conclude that two xabs components are suf- 
ficient to model the warm absorption in ESO 1 13-GOlO. Finally, 
we repeated the above analysis using other SEDs adopted for 
ESO 1 13-GOlO, which are described in Sec.|5]and shown in Fig. 

m 

The best-fit parameters of the model (obtained using SED 
E) are shown in Table [1] Note that the velocities of the warm 
absorber phases are the same for the different SEDs of Fig. |4] 
Furthermore, the only parameter of the warm absorber outflow s 
which significantly changes as a consequence of using different 
SEDs is the ionisation parameter ^, whose values are given in 
Fig. |5] The total column densities of different ions in the two 
warm absorber phases are listed in Table |2] In the top panel of 
Fig. |2] the RGS spectrum and the best-fit model are presented 
and in the bottom panel we show how each phase in our model 
contributes to the absorption. In Fig.|3] we show close-ups of the 
RGS spectrum where the most prominent absorptions lines are 
detected. 



5. Modelling with different Spectral Energy 
Distributions 

In order to calculate a more accurate ionisation balance required 
to improve the photoionisation modelling of the warm absorber, 
it is essential to determine the broad-band continuum of the 
source. Note that there are few flux measurements for ESO 113- 
GOIO compared to some of the well-studied Seyfert galaxies. In 
the following subsections we describe how different parts of the 
SED were constructed. In Sect. I5.4l we discuss the uncertainties 
in the construction of the SED due to nuclear obscuration of the 



Table 1. Best-fit parameters of the power-law continuum and 
two-phase warm absorber xabs model obtained from fitting the 
RGS spectrum as described in Sect. |4]and using SED E in Fig. 
m The;^J = 1.12 (1046 d.o.f.). The fitted parameter errors are 
quoted at 90% confidence for one interesting parameter 



Power-law: 


2.62 + 0.05 






Normj^ 
3.1+0.1 




xabs 


log^^ A'h 


d 


Flow 


vl RMS 




Phase 1: 


3.2 + 0.1 1.6 ± 


1.0 


-1100 = 


t 300 500 ± 


300 


Phase 2: 


2.3+0.1 0.5 ± 


0.2 


-700 ± 


100 100 ± 


40 



" Photon index. 

* Normalisation in 10" photons s"' keV ' at 1 keV. 

erg cm s"'. 

10^2 cm-^ 
' kms"'. 



source and how, in order to investigate the effects of these uncer- 
tainties, we select different SEDs in our further modelling. 

5.1. Interstellar de-reddening of our Galaxy 

We corrected all the optical-UV fluxes for interstellar reddening 
in our Galaxy using the reddening curve o f^ardelli etal.(198j), 
including the update for near-UV given by lO'DonnelU (11994b . 
The Galactic interstellar c olour ex cess E (B - V ) = 0.025 mag 
is based on calculations of iSchlegeletal.! (119981) as shown in the 
NASA/IPAC Extragalactic Database (NED). Also in this case Ry 
was fixed at 3.1. 

5.2. Optical-IR-Radio part of the SED 

Optical flux measurements in Cousin's B (440 nm) and R 
(640 nm) bands were taken from the Surface Photometry 
Catalogue of the ESO-Uppsala Galaxies dLauberts & ValentiinI 
1989i) . We extracted IR flux measurements at J (1.25 //m), H 
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Table 2. The total column densities of the most relevant ions 
of the warm absorber outflow, derived from the xabs modelling 
described in Sect. |4] Percentages of ionic column density pro- 
duced by each warm absorber phase are also shown. The ranges 
of values given correspond to the minimum and maximum val- 
ues found using the different SEDs of Fig.|4] 
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(1.65 /Km) and Ks (2.17 u rn) bands from fli e 2MASS Afl-Sky 
Catalog of Point Sources (ICutri et al.l l2003h . Furthermore, we 
used far-IR data at 12 fim, 25 yum, 60 i^m and 100 //m from the 
IRAS Faint Source Catalog (Moshiret al. 1990), and at 160 yum 
from die AKARI/FIS All-Sky Survey Bright Source Catalog 
(lYamamura et al.ll2010h . These data points are shown as circles 
in the SEDs displayed in Fig. |4l the red circles are for fluxes 
only corrected for Galactic extinction as described in Sect. 15.11 
whereas the fluxes shown as blue circles are also corrected for 
extinction by the host galaxy of the AGN as described in Sects. 
[37I](SEDs C and D) and [32] (SEDs E and F). We note that the 
fluxes used here are AGN-dominated but include some stellar 
optical emission coming from the nuclear region of the galaxy 
and IR dust emission from the whole galaxy; however, for the 
purpose of our study, contaminating emission by the host galaxy 
would have a negligible effect on the broad-band SED and pho- 
toionisation modelling. 

The only radio-band flux measurement available for 
ESO 113-GOlO was found in the Sydney University Molonglo 
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Fig. 2. Top panel: RGS spectrum of ESO 1 13-GOlO, fitted using 
the two-phase xabs model described in Sect.|4|and SED E in Fig. 
m The data are shown in red and the model in blue. Residuals 
of the fit, (Observed-Model)/Model, are displayed below the 
panel. Bottom panel: The absorption components of the best-fit 
RGS model as described in Sect. |4] with the best-fit parame- 
ters of the model given in Table [1] The blue component shows 
the Galactic absorption applied to the power-law continuum (in 
black). The two warm absorber phases (Phase 1 in red and Phase 
2 in green) are applied separately to the Galactic-absorbed con- 
tinuum, showing how each phase contributes to the absorption. 

Sky Survey (SUMSS) Source Catalog (iMauch et al.ll2003h . The 
flux at 36 cm (outside the plots in Fig.|4]l is 13.9 + 1 .0 mJy. 

5.3. UV-X-ray pan oi the SED 

The UV part of the SED was constructed from the XMM-Newton 
OM observation and the GALEX All-sky Imaging Survey (AIS). 
The 2005 OM data were taken in the U (3440 A), UVWl 
(2910 A), UVM2 (2310 A) and UVW2 (2120 A) filters. The 
2007 GALEX data were taken in the NUV (2267 A) and FUV 
(1516 A) filters. Note that the OM UV data (which ai-e simulta- 
neous with the X-rays) are not simultaneous with the GALEX 
UV data. However, the OM UVM2 flux is consistent with the 
GALEX NUV flux, which indicates that the UV flux was at 
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Fig. 3. Close-ups of parts of the RGS spectram and the best- 
fit model (the two-phase xabs modelling described in Sect. |4] 
and using SED E in Fig. |4]i where the most prominent absorp- 
tion lines are detected. The data are shown in red and the model 
in blue. Residuals of the fit, (Observed-Model)/Model, are dis- 
played below each panel. 



a similar level during the XMM-Newton and GALEX observa- 
tions. The OM and GALEX data are shown as squares in the 
SEDs displayed in Fig.HJ the red squares show fluxes only cor- 
rected for Galactic extinction as described in Sect. 15.11 whereas 
the blue squares are for fluxes also corrected for extinction in the 
host galaxy of the AGN as described in Sects. [3T|(SEDs C and 
D) and[l2](SEDs E and F). 

The 0.3-10 keV X-ray part of the SED was constructed 
using the EPIC-pn data, corrected for the Galactic Hi col- 
umn den sity in our line of s ight, A ^h = 2.78 x 10^" cm"^, 
given by iDickev & LockmanI (Il990h . The EPIC-pn data are 
shown as small green circles in the SEDs displayed in Fig. |4] 
Unfortunately, there are no X-ray data available for ESO 113- 
GOIO above 10 keV; therefore in order to estimate the SED 
above 10 keV, we extrapolated the high energy part of the EPIC- 
pn spectrum using a power-law model up to 20 keV (4.84 x 10'^ 
Hz); then from 20 keV to 162 keV, th e INTEGRAL data from 
the SED of Mrk 509 as given in lKaastra et al. (201 1.) were used, 
but scaled down to match the ESO 1 1 3-GO 1 flux at 20 ke V. This 
way the shape of the ESO 1 13-GOlO SED above 10 keV resem- 
bles the Compton 'hump' commonly seen in Seyfert AGN. We 
have tested whether the shape of the continuum above 10 keV, 
and in particular the shape and flux of the Compton hump, af- 
fect the results of the warm absorber analysis and the thermal 
stability curves (Sect. 15.5b : we find that changing the hard X- 
ray continuum does not modify our results and all the derived 
parameters remain unchanged within the given errors. 

5.4. Different assumed SEDs for modelling 

Different SEDs of the ionising source can have a significant ef- 
fect on the ionisation balance of the absorbing material and con- 
sequently aff'ect the stability and structure of the warm absorber 
So we have taken into account the uncertainties by adopting dif- 
ferent SEDs in our further modelling of the warm absorber 

For the case of ESO 1 13-GOlO, the largest uncertainty arises 
from the level of extinction and reddening in the host galaxy of 
the AGN, reported in Sect.|3] (1) there is uncertainty in measur- 
ing the Balmer decrement in ESO 1 13-GOlO and (2) the Balmer 
decrement can be intrinsic to the BLR or NLR clouds, so there 
are uncertainties on deriving reddening information from the 
Balmer decrement. As shown in Fig.|4] we considered SEDs cor- 
rected only for Galactic extinction (SEDs A, B), and with also 
the AGN obscuration taken into account as described in Sect. 
l3.1l (SEDs C, D) based on the Balmer decrement and Sect. 13.21 
(SEDs E, F) based on the continuum reddening. Furthermore, 
the precise locations of the warm absorber outflows relative to 
the putative dusty torus are unknown. A dusty torus is expected 
to exist around the nuclear source in AGN and be responsible 
for the IR 'bump' seen in their SEDs by reprocessing of the nu- 
clear high energy emission. So there is uncertainty as to how 
much I R radiation is s een by the warm absorber outflows. In ad- 
dition, Poiwietitin (|20075 report that the 2001 EPIC-pn data 
show no significant constant narrow 6.4 keV Fe Ka line (with 
an upper limit of EW < 32 eV), hence suggesting lack of any 
dominant Compton reflection emission from distant cold matter 
such as the dusty torus. Thus, to investigate the eff'ect of emis- 
sion from the dusty torus in our ionisation balance calculations, 
we consider the same SEDs as described above but without the 
IR bump (SEDs B, D, F, H in Fig.|4| to represent cases in which 
the outflow does not see emission from the dusty torus. 

We adopted these SEDs in our calculations of the ionisation 
balance and thermal stability curves, as described in Sect. 15.51 
We note here that adopting different SEDs does not affect the 
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Fig. 4. The different SEDs (A-H) assumed in the modelling of ESO 113-GOlO warm absorber as described in Sects. |4]and|5] which 
were used in the ionisation balance calculations. The IR and optical data points (Sect. I5.2l i are shown as circles and the UV data 
(Sect. 15.3b as squares. The data points in red are corrected only for Galactic extinction as described in Sect. 15.11 whereas the blue 
data points (in panels C, D, E and F) are also corrected for nuclear obscuration in the AGN as described in Sects. l3.1l and U!2l The 
SEDs in panels C and D have been corrected for nuclear obscuration by using the Balmer decrement and the ones in panels E and 
F have been corrected for nuclear obscuration by using the continuum reddening. In panels C and D, the SEDs with some parts 
shown in grey and no data points correspond to nuclear obscuration correction with an intrinsic E{B -Y)- 0.64 (named SEDs 
CI and Dl) and the ones with blue data points superimposed on the curve correspond to nuclear obscuration correction with an 
intrinsic E{B -Y) - 0.73 (named SEDs C2 and D2). The radio point (Sect. 15.2b is not displayed in this figure due to its low flux. 
The EPIC-pn X-ray data corrected for Galactic absorption are shown as small green circles. The black diamonds correspond to flux 
at 20 keV extrapolated from the EPIC-pn spectrum, and the extrapolations to the higher energy parts of the SED, shown as dotted 
lines, are described in Sect. 15.31 SED G is the SED of Mrk 509 introduced in Sect. |4] and used as a 'standard' representation for a 
Seyfert AGN. SEDs on the right-hand side (B, D, F, H) are the same as their adjacent SEDs on the left (A, C, E, G) but without the 
IR bump and are used in modelUng cases where the warm absorber does not receive IR radiation from the dusty torus. 
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goodness of the RGS fits since the columns of the relevant ions 
detected by the RGS are not significantly different for different 
SEDs. The parameter which changes significantly as a result of 
using different SEDs is the ionisation parameter ^, which is given 
in Sect. l5.5l (Fig.l5]) for each case and plays a significant role in 
determining the structure of the warm absorber using thermal 
stability curves. 

5.5. Different thermal stability curves for the warm absorber 

The ionisation balance required for the photoionisation mod- 
elling is determined by the SED. We now study the structure of 
the warm absorber for each corresponding SED case in Fig. 
To investigate the stability of each warm absorber phase and the 
possibility of the two phases co-existing in pressure equilibrium, 
the pressure form of the ionisation parameter, H (introduced by 
iKrolik et aLlfTosTh . needs to be used. The parameter H is defined 
as 



(4) 



Ancr^nkT 



where c is the speed of light, k the Boltzmann constant and T the 
gas temperature (and L, r and n as in Eq.|3]l. Substituting Eq. ^ 
into Eq. (|4]i gives 



AnckT 



19222 - 
T 



(5) 



For each SED we produced the corresponding thermal sta- 
bility curve, shown in Fig.|5] using the output of the Cloudy runs 
(see Sect. |4]i. The best-fit ionisation parameters of Phase 1 and 
Phase 2 and their errors are marked on the stability curves. In 
order for the two phases to co-exist in pressure equilibrium, they 
must have overlapping values of H. We discuss these results in 
Sect.lO 



6. Discussion 

6.1. Nuclear obscuration or an intrinsically large Balmer 
decrement? 

Here we discuss several possibilities for the discrepancy between 
the large Balmer decrement in ESO 1 13-GOlO and the fact that 
the X-ray spectru m is not intrinsically absorbed by cold neutral 
gas. We note that i Porquet et al.l (12007 ) obtained an upper-limit 
to the neutral gas column density in ESO 1 13-GOlO from anal- 
ysis of the EPIC-pn spectrum, which is much lower than that 
predicted for a typical Seyfert 1.8 galaxy, in which the BLR is 
expected to be reddened. Porquet et al. (2007) briefly mentioned 
possible explanations for such a discrepancy (explanations 1, 2 
and 4 which we expand on below), but did not go into details. In 
this work, with the benefit of our simultaneous optical/UV (OM) 
and X-ray (RGS, EPIC-pn) analysis of the XMM-Newton data 
and also analysis of the archival optical spectrum of ESO 113- 
GOIO, we can discuss the likelihood of each explanation and sug- 
gest which is the most viable solution. 



(1) Non-simultaneous optical and X-ray/UV observations 

The optical spectrum in which ESO 113-GOlO shows the large 
Balmer decrement was taken in 1996, whereas the XMM-Newton 
observations were taken in 2001 and 2005. The nearest X- 
ray data to the time of the optical observation were obtained 
in 1995 by ROSAT. So we can compare the X-ray flux in 



2001 and 2005 with the flu x in 1995 to check flie possibil- 
ity of variability. As given in iPietsch et al.l (Il998l) . the soft X- 
ray (0.1-2.4 keV) energy flux from the 1995 ROSAT obser- 
vation is 8.5 X 10"'^ erg s"' cm"^. We calculate the soft X-ray 
flux in the same band to be 3.7 x 10"'^ erg s"' cm"^ in 2001 and 
5.8 X 10^'^ erg s"' cm"^ in 2005 from the XMM-Newton obser- 
vations. This flux change over the years is too small to be at- 
tributed to occultation by the column of gas inferred from the 
Balmer decrement in Sect. |33] A^h =(3.5-5.0) xlO^' cm"^ 
Assuming such a column density of cold gas was present dur- 
ing the optical observation and the unabsorbed soft X-ray flux 
was at the same level as in 2001 or 2005 (when we know there 
was no intrinsic neutral absorption) would imply an obscured 
soft X-ray flux smaller by a factor of 3.1^.1, whereas in 1995 
the X-ray flux was actually larger by a factor of 1 .5-2.3. 

More importantly, during the 2005 XMM-Newton observa- 
tion, when simultaneous UV (4 UV filters of OM) and X-ray 
observations are available, the optical-UV continuum, which is 
supposed to represent thermal emission from the accretion disc, 
appears to be intrinsically reddened as described in Sect. 13.21 So 
the non-simultaneity of the optical spectroscopy and X-ray ob- 
servations is unlikely to be the cause of discrepancy between the 
reddening and X-ray neutral absorption. 



(2) Partially covering neutral gas clouds 

The clouds of neutral gas and dust that obscure the NLR and 
BLR may only partially cover the nuclear source, such that they 
are not in our line of sight to the X-ray emitting region. This way 
optical emission from the NLR and BLR is reddened but the X- 
rays remain unabsorbed. This scenario is however an unlikely 
possibility as it requires an ad hoc geometry of the absorbing 
regions. Furthermore, since the slope of the optical/UV contin- 
uum shows that intrinsic UV emission from the accretion disc 
has been reddened, the intrinsic X-ray emission originating from 
the corona in close vicinity of the disc must also have passed 
through the same material as the UV emission. Therefore, par- 
tially covering neutral gas clouds in the host galaxy of the AGN 
are not a feasible solution. 



(3) Intrinsically large Balmer decrements 

One possible explanation is that Balmer decrements are intrin- 
sically large for the BLR or even the NLR. The physical condi- 
tions in the BLR are to some extent uncertain and higher density 
clouds within the BLR can cause collisional and radiative effects 
which can rai se the value of the int rinsic Balmer decrement, e.g. 
from 3 to 10 (iKwan & KroUklfigSli) . 

iBarcons et al.l ( l2003h have analysed the optical (WHT) and 
X-ray {XMM-Newton EPIC-pn) spectra of the Seyfert 1.8/1.9 
galaxy HI 320+55 1 . They find that the Balmer decrement Ha/HyS 
is about 6 for the NLR and about 27 or more for the BLR. 
Despite such large Balmer decrements, the X-ray spectrum 
is not intrinsically absorbed. They find a 3cr upper limit of 
1.4 X 10^*^ cm"^ to the obscuring column density, which they 
report is about 70 and 7 times smaller than the minimum pre- 
dict ed from th e BLR and NLR Balmer decrements respectively. 
Bar cons et al. I (l200l also rule out the existence of an ionised 
absorber in H1320+551 from their X-ray spectral fits. They sug- 
gest the large Balmer decrement of the BLR is an intrinsic prop- 
erty and not caused by internal reddening, and so H1320-I-551 is 
not consistent with being an obscured Type-1 Seyfert AGN. 
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Fig. 5. Different stability curves (A-H) calculated for the corresponding SEDs (A-H) of Fig.H] In panels C and D, the curves in 
blue correspond to SEDs CI and Dl and the ones in black to SEDs C2 and D2. The positions of the warm absorber phases (Phase 1 
and 2), along with their 2cr (95.4%) errors, are indicated as thick magenta strips on the curves. Phase 1 is the one with higher T . The 
ionisation parameters of Phase 1 and 2 (i.e. ^\ and ^2) in units of erg cm s"' are given for each corresponding stability curve. The 
dotted lines are plotted for display purposes to show the range of H for one of the phases on each curve. The regions overlapping in 
H on each stability curve are in pressure equilibrium. 



However, in ESO 113-GOlO, the UV continuum provides 
an independent estimate of the reddening from that inferred 
from the Balmer decrement. As shown in Sect. 13.2] the amount 
of reddening derived from the continuum study, E{B - V) - 



0.39, is less than that inferred from the Balmer decrement, 
E(B-V) - 0.64-0.73. This may indicate that the Balmer decre- 
ment over-estimates the amount of reddening, as the intrinsic 
Ha/HyS is larger than the theoretical values 3.1 -3.5. From Eq. 
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|2] the intrinsic Ha/H/3 needs to be about 4.8 in order to imply 
E{B - y) = 0.39, while the value observed is ~ 8. One rea- 
son disfavouring the large reddening obtained from the Balmer 
decrement is that the bolometric luminosity calculated from the 
de-reddened SEDs (SEDs CI and C2 in Fig. lU is high for a 
Seyfert galaxy: 2.8-5.8 xlO^^ erg s"'; on the other hand, the 
bolometric luminosity of the SED de-reddened using the con- 
tinuum (SED E in Fig.|4]i is 4.7 x 10'*^ erg s"' . This is very simi- 
lar to the bolometric luminosity of Mrk 509, calculated from the 
SED G: 4.6 x 10"^^ erg s"' . Nonetheless, even if to some degree 
the large Balmer decrement is intrinsic, the slope of the optical- 
UV continuum indicates there must be intrinsic reddening taking 
place in our line of sight in ESO 113-GOlO. 

(4) Dusty warm absorber 

The most viable explanation for the discrepancy between the 
reddening and X-ray absorption in ESO 113-GOlO is that the 
X-ray warm absorber contains dust. The dust causes the ob- 
served optical-UV reddening of the continuum and the large 
Balmer decrement, whilst the X-rays remain photoelectrically 
unabsorbed since the gas in the warm absorber is ionised (pre- 
dominantly due to the ionisation of carbon and oxygen), and the 
X-ray opacity of dust is lower than that of cold neutral gas. This 
means unusually large dust-to-gas ratio clouds in the ISM of 
the host galaxy are not needed to explain the discrepancy be- 
tween reddening and X-ray absorption. Previous studies show 
that dust particles can survive the dust destruction mechanisms 
(sublimation and thermal sputtering) under conditions fo und in 
some warm absorbers. For example, Revn olds et al.l(ll997 ) sug- 
gest the existence of a dusty warm absorber in the Seyfert- 1 
galaxy MCG-6-30-15, to account for the discrepancy between 
the reddening and lack of X-ray absorption by neutral gas, and 
discuss the survival of dust grains under conditions found in a 
photoionised warm absorber. We check here if dust can survive 
in the warm abso rber of ESQ 1 1 3- GO 10. 

As shown in iBarvainisI (Il987h . the sublimation radius (the 
minimum distance from the central source at which particular 
grains can exist) for graphite grains is given by 

^sub = 1.3^46 71-5^0^ PC (6) 

where Luv,46 is the UV luminosity of the central source in units of 
10"*^ erg s"' , and Tisoo is the graphite grain sublimation tempera- 
ture in units of 1500 K. Taking an upper estimate of lO'** erg s"' 
for the UV luminosity of ESO 113-GOlO (consistent with our 
SEDs), we obtain /?sub ~ 1.3 pc. So as long as the dust resides 
in a warm absorber phase which is further than 1 .3 pc from the 
source, it is not destroyed by sublimation. 

Another mechanism that can destroy dust grains is sput- 
tering, in which atoms or molecules are knocked off the sur- 
fa ce of th e dust grain due to collisions with hot ions. As given 
in iBurke & Silk (1974), the threshold for sputtering graphite 
grai ns corresponds to a gas temperature of about 4 x 10^ K. 
iBurke & Silk (1974J have shown that above the sputtering 
threshold, the grain lifetime against destruction by sputtering is 

f,pu, = 6.25xlO>'(Fnr°-^)-'s, (7) 

where Y is the sputtering yield at T4, the incident particle en- 
ergy in units of 10"* K, and n is the incident particle density. 
Phase 2 of the warm absorber in ESO 113-GOlO has a max- 
imum temperature of 2.6 x 10^ K corresponding to SED CI 
(see Fig. |5]l. This temperature is below the sputtering thresh- 
old and thus graphite dust grains in Phase 2 should not be de- 
stroyed by the sputtering mechanism. On the other hand. Phase 



1, which is hotter than Phase 2, has a minimum temperature 
of 9.3 x 10^ K (corresponding to SED G) and a maximum tem- 
perature of 2.7 X lO*" K (corresponding to SED A). So the tem- 
perature of Phase 1 is higher than the sputtering threshold for 
all SED cases. Substituting n from the expression of the ioni- 
sation parameter (^ - L/nr^) into that of fsput (Eq. |7]), we find 
that for Phase 1 (with the parameter values given in Table [1}, 
fsput ~ 22.8 yr, where r is the distance of the gas from the cen- 
tral source in pc. Furthermore, we can roughly estimate Phase 1 
flow timescale, fflo„ ~ r/v, where v is the velocity of Phase 1 (i.e. 
1 100 km s~'); so fflow ~ 890 r yr, where the distance r is in pc. 
The dust grains are destroyed by sputtering if fsput < fflow; there- 
fore, if r < 39 pc, then the dust grains in Phase 1 cannot survive 
the sputtering. Assuming a constant velocity outflow, an upper 
limit for the distance r can be estimated by using A^h = nrf 
to substitute into ^ - L/nr^ and eliminate «, where we take the 
volume filling factor / = 1 . For Phase 1 we obtain r < 52 pc, so 
the dust grains can only survive in Phase 1 if 39 < r < 52 pc. 
However, Phase 2 (with r < 1.3 kpc applying the same calcu- 
lations as above), which is the lower-ionised phase of the warm 
absorber with a temperature below the sputtering threshold, is 
more likely to host dust grains than Phase 1 . 

6.2. The structure of the warm absorber 

From modelling of the RGS spectrum, we found that the warm 
absorber in ESO 113-GOlO consists of at least two phases of 
ionisation. As shown in Table [1] the higher-ionised Phase 1 has 
a larger column density than the lower-ionised Phase 2, and also 
has larger outflow and RMS velocities. The thermal stability 
curves of Fig. |5] show that depending on which ionising SED is 
selected, quite different stability curves are obtained. The SED 
also dictates whether the absorber phases lie on the stable or un- 
stable part of the curves. Only in case C, in which the optical-UV 
continuum is de-reddened using the Balmer decrement and the 
warm absorber does receive IR radiation from the dusty torus of 
the AGN, the two phases have overlapping values of S and hence 
can be in pressure equilibrium. However, for all other cases, the 
two phases do not overlap in H, so they are not in pressure equi- 
librium. Since also their velocities are different, we conclude that 
the two phases are likely to be separated and out of pressure 
equilibrium. 

The exact location of the warm absorber relative to the dusty 
torus is unknown, so it is not clear whether or not the warm ab- 
sorber sees IR radiation from the toru s. As mention e d in S ect. 
15.41 from analysis of the Fe Ka line, iPorquet et al.l (l2007h re- 
port lack of Compton reflection from distant cold regions such 
as the dusty torus. However, lack of Compton reflection from 
cold matter does not imply absence of a dusty torus, because the 
BLR is likely to be obscured by the torus given the lack of broad 
optical emission lines in this object; also it does not imply the 
warm absorber is not receiving IR emission from the torus be- 
cause the warm absorber phases (most likely evaporated off the 
torus itself), are likely to be further away from the torus (see the 
r values derived in the previous section), so they may see IR ra- 
diation from the outer parts of the torus. As inferred from Fig.|5] 
apart from cases C and D, the IR radiation from the dusty torus 
does not affect whether or not the two phases are in pressure 
equilibrium anyway. 

As discussed earlier, dust can survive in the lower-ionised 
phase (Phase 2) of the warm absorber, which can account for 
the discrepancy between the large optical/UV reddening and the 
low amount of X-ray neutral absorption. Since the NLR is red- 
dened as inferred from the Balmer decrement, then the dust in 
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Phase 2 of the warm absorber cannot be closer to the nucleus 
than the NLR. So Phase 2 of the warm absorber must lie some- 
where between the NLR and the maximum distance of 1 .3 kpc at 
which the absorber can be kept photoionised. Spatially resolved 
NLR in nearby S eyferts show diameters of order 10^-10^ pc 
dOsterbrock & Fer land 2006). The range of possible distances 
from the nucleus for the dusty warm absorber (Phase 2) overlaps 
the range which the NLR is likely to occupy, so the two regions 
could be co-spatial. On the other hand, the higher-ionised Phase 
1, which is not co-spatial with Phase 2 on the grounds of stabil- 
ity analysis and different outflow velocities, has r < 52, which 
places it interior to Phase 2 and the NLR. 



7. Conclusions 

We have presented the first-ever study of the X-ray warm ab- 
sorber and nuclear obscuration in the Seyfert 1.8 ESO 113- 
GOIO. We conclude fliat: 

1 . The warm absorber which is detected by the XMM-Newton 
RGS consists of two phases of ionisation, with log^ ^3.2 
and 2.3 respectively. The higher-ionised component (Phase 
1), with A^H ~ 1.6 X 10^^ cm"^, is outflowing with v w 
-1 100 km s '; the lower-ionised component (Phase 2), with 
A^H ~ 0.5 X 10^- cm"^, is outflowing with v w -700 km s"'. 

2. ESO 113-GOlO displays a large Balmer decrement 
(Ha/HyS ~ 8) relative to the theoretical predictions for the 
NLR and BLR, implying a significant amount of reddening 
with E{B -V)- 0.64-0.73 and a large column of cold gas, 
with A^H =(3.5-5.0) xlO^' cm"^, in the host galaxy of the 
AGN. However, the X-ray spectrum is not absorbed by such 
column of gas; we find an upper limit of 8.9 x 10'^ cm"^ for 
the neutral hydrogen equivalent column density A^h- 

3. Alternative explanations for the discrepancy between the 
large Balmer decrement and lack of X-ray absorption are dis- 
cussed. We find partially covering gas clouds or an intrinsi- 
cally large Balmer decrement are not feasible solutions since 
the shape of the optical-UV continuum flux (supposed to be 
thermal emission from the accretion disc) also indicates red- 
dening. We suggest a dusty warm absorber is a viable expla- 
nation for the discrepancy, without requiring unusually large 
dust-to-gas ratio clouds in the ISM of the host galaxy. 

4. We show that graphite dust grains can survive the sublima- 
tion and sputtering processes in the lower-ionised phase of 
the warm absorber (Phase 2), and thus can cause the ob- 
served optical-UV reddening, whereas the X-rays remain un- 
absorbed due to lack of neutral column of gas in the ionised 
absorber 

5. We have explored the uncertainties in the construction of the 
SED of ESO 113-GOlO due to nuclear obscuration of the 
source and the effects these have on the results of our warm 
absorber analysis and the ionisation balance calculations re- 
quired for photoionisation modelling. The only parameter of 
the warm absorber outflows which significantly changes as 
a consequence of using different SEDs is the ionisation pa- 
rameter Furthermore, different SEDs result in different ther- 
mal stability curves, which determine whether or not the two 
phases of the warm absorber are in pressure equilibrium. We 
find that only for the case where the SED is corrected for 
reddening using the Balmer decrement in the host galaxy and 
the warm absorber does receive IR radiation from the dusty 
torus expected to surround the AGN, the two phases of ion- 
isation can co-exist in pressure equilibrium on the thermal 



stability curve, whereas for other SEDs the two phases are 
distinct and not in pressure equilibrium. 
6. This work demonstrates the importance of establishing the 
intrinsic broad-band continuum SED in order to interpret 
correctly the structure of the warm absorber We conclude 
that the two phases of the warm absorber in ESO 1 13-GOlO 
are most likely separate distinct phases and not in pressure 
equilibrium; and the dusty phase of the warm absorber is 
likely to be co-spatial with the NLR. 
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